Properties of the multicritical point of ±J Ising spin glasses on the square lattice 
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We use numerical transfer-matrix methods to investigate properties of the multicritical point of 
binary Ising spin glasses on a square lattice, whose location we assume to be given exactly by a 
conjecture advanced by Nishimori and Nemoto. We calculate the two largest Lyapunov exponents, 
as well as linear and non-linear zero-field uniform susceptibilities, on strip of widths 4 < L < 16 
sites, from which we estimate the conformal anomaly c, the decay-of-correlations exponent n, and 
the linear and non-linear susceptibility exponents j/v and j nl /v, with the help of finite-size scaling 
and conformal invariance concepts. Our results are: c = 0.46(1); 0.187 < r\ < 0.196; y/u = 
1.797(5); 7™ /v — 5.59(2). A direct evaluation of correlation functions on the strip geometry, and 
of the statistics of the zeroth moment of the associated probability distribution, gives r\ = 0.194(1), 
consistent with the calculation via Lyapunov exponents. Overall, these values tend to be inconsistent 
with the universality class of percolation, though by small amounts. The scaling relation y nl jv — 
2-y/u + d (with space dimensionality d = 2) is obeyed to rather good accuracy, thus showing no 
evidence of multiscaling behavior of the susceptibilities. 

PACS numbers: 75.50.Lk,05.50.+q 



I. INTRODUCTION 

The critical behavior of Ising spin glasses has been the 
subject of intensive investigation in the recent pasfci. A 
number of results have been derived, both analytically 
and numerically; however, many aspects of interest have 
not been fully elucidated so far. Among these are the 
properties of the multicritical point which is known to 
exist for suitably low concentrations of antiferromagnetic 
interactions, even for two-dimensional lattices in which 
a spin glass phase is not expected to occur at non-zero 
temperatures. 

Here we consider binary (±J) spin glasses, i.e., we 
assume that Ising spin- 1/2 magnetic moments interact 
via nearest-neighbor couplings Jjj of equal strength, and 
whose signs are given by the quenched probability distri- 
bution: 

P(J l3 ) = p5(Jij - Jo) + (1 - p) 5(Jij + Jo) . (1) 

In this case, the phase diagram on the temperature- 
concentration (T — p) plane exhibits a critical line which, 
for low concentrations 1 — p of antiferromagnetic bonds, 
separates ferro- and paramagnetic phases-- 2 . As p de- 
creases, so does the transition temperature. Below a 
critical concentration p c , ferromagnetic order disappears, 
and a spin glass phase emerges at T = 0. For space di- 
mensionality d > 3, the spin glass phase extends to finite 
temperatures as well. The Nishimori line (NL) is a spe- 
cial line on the T — p plane, defined by: 



-2J a /T 



1-p 



(NL, p> 1 -) 



(2) 



On this line, several exact results have been obtained^. In 
particular, the configurationally averaged internal energy 
is an analytical function of T, even at the multicritical 
point (the Nishimori point, NP) where the NL crosses the 
ferro-paramagnetic phase boundary^. Furthermore, the 
NL is invariant under renormalization-group transforma- 
tions, so the NP corresponds to a fixed point. Numerical 
work in two-dimensional systems^ shows that internal- 
energy fluctuations along the NL go through a maxi- 
mum at the NP, thus indicating that the latter indeed 
marks a change in the behavior of the distribution of 
frustrated plaquettes. This, in turn, is consistent with 
the picture that the phase transition at the NP is of 
geometry-induced nature^ (though not necessarily in the 
same universality class of random percolation) . 

Recently it was predicted^ that, on a square lattice, 
the NP should belong to a subspace of the T — p plane 
which is invariant under certain duality transformations. 
For ± J Ising systems the invariant subspace is given by^ 

p\og 2 (l + e- 2J o/ T ) + (l-p) log 2 (l + e 2J °/ T ) = i . (3) 

The intersection of Eqs. <|2jl and (pj gives the conjec- 
tured exact location of the NP, namely p = 0.889972 • • • , 
T/J = 0.956729 • • • , to be referred to as CNP. In pre- 
vious work2A&, approximate estimates for the location 
of the NP were used in the calculation of the associated 
critical exponents, with the overall conclusion that the 
transition there does not belong to the universality class 
of random percolation. A numerical study of correlation- 
function statistics at the CNR& points to a similar pic- 
ture. 
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Extensions of the conjecture to triangular and honey- 
comb lattices have been proposed 11 ^, and verified by nu- 
merical studies to a fairly good degree of accuracy^. Ev- 
idence thus far available indicates that the critical prop- 
erties of the NP in two-dimensional ±J Ising systems 
are universal in the expected (i.e., lattice-independent) 
sensed, though they belong to a distinct universality 
class from that of percolation. 

Here we use numerical transfer-matrix methods, to- 
gether with finite-size scaling and conformal invariance 
concepts, to investigate critical properties of the NP of 
± J Ising spin glasses, on long strips of a square lattice. 
We shall assume the location of the multicritical point to 
be that of the CNP given above. Indeed, previous work 
(see Ref.^and references therein) strongly indicates that, 
even though the conjecture may turn out not to be exact, 
it is certainly a very good approximation to the actual 
position of the NP. 

In Section [H] we evaluate the central charge, or con- 
formal anomaly^ 3 -. As this is given by the coefficient of 
the finite-size correction to a bulk quantity (the criti- 
cal free energy) which is itself not known exactly for 
the present case, one has many sources of uncertainty 
to contend with, not to mention those intrinsic to the 
sampling of quenched disorder configurations. By work- 
ing at the CNP, we attempt to eliminate one such source 
which is the location of the critical point. In the absence 
of further exact results, whether or not such choice in 
fact introduces systematic distortions can only be found 
by comparative analysis of numerical data pertaining 
to the problem. In Section IIIII we calculate the decay- 
of-correlations exponent related to the zeroth moment 
of the correlation-function probability distribution, both 
via the difference between the two largest Lyapunov ex- 
ponents, and by direct evaluation of correlation func- 
tions as done in earlier work^iS. This, together with 
the use of pertinent conformal-invariance relationships, 
yields further independent evidence related to the uni- 
versality properties of correlations at the NP. In Sec- 
tional both linear and non-linear zero-field susceptibil- 
ities are investigated. While the former have been eval- 
uated previously for square 7 -, as well as triangular and 
honeycomb^, lattices, no results for the latter appear 
to be available. As explained below, the scaling of non- 
linear susceptibilities may give indications of multiscaling 
behavior. Finally in Section EI concluding remarks are 
made. 



II. FREE ENERGY AND CENTRAL CHARGE 

We consider strips of width L sites and periodic bound- 
ary conditions across. For consistency with earlier work 7 , 
we used only even widths, in order to accommodate pos- 
sibly occurring unfrustrated antiferromagnetic ground 
states (though later results showed in practice that, at 
least for the relatively low concentrations of antiferro- 
magnetic bonds around the NP, no noticeable distortions 



arise when odd values of L are considered as well2ii£). 
Appropriate sampling of quenched disorder is produced 
by using strip lengths N 3> 1, along which bond config- 
urations are drawn from the distribution Eq. QJ. The 
configurationally-averaged (negative) free energy /z, (in 
units of ksT) is given by 



f L =L- l K {L) 



(4) 



where Ao is the largest Lyapunov characteristic expo- 
nent, extracted from the product of N — » oo transfer 
matrices Tj which connect site columns j and j + i.e., 



A = lim — In 



N 



(5) 



where |i>o) is an arbitrary initial vector of unit modu- 
lus. Higher-order exponents may be obtained through 
iteration of a set of initial vectors orthogonal both 
mutually and to \vq), with adequate reorthogonalization 
every few steps, to avoid contamination^. 

In our calculations we have used N ~ 10 5 — 10 6 . 
The uncertainty related to the finite number of terms 
in Eq. © is estimated as follows. In order to avoid tran- 
sient effects, the first Ao ~ 10 3 iterations are discarded. 
Accumulated averages are evaluated and stored, for each 
10 3 subsequent iterations. From this set of averages, one 
then calculates global averages and their corresponding 
fluctuations. In our study, we have always made use of 
a canonical distribution of disorder, that is, the + Jo and 
— Jo couplings are randomly extracted from a reservoir 
which initially contains exactly as many of each as given 
by Eq. with the value of p corresponding, e.g., to 
the CNP. This way, sample-to-sample fluctuations are 
considerably reduced. Final estimates of the free energy 
and other quantities of interest have been extracted from 
arithmetic averages for distinct disorder realizations. 

The conformal anomaly, or central charge c, which 
characterizes the universality class of a conformally- 
invariant model at the critical point, can be evaluated 
via the finite-size scaling of the free energy on a strip 
with periodic boundary conditions across^, 



f(T c ,L) = f(T c ,oo) 



7TC 

6L2 



O 



(6) 



where /(T c ,oo) = lirn^oo f(T c ,L) is a regular term 
which corresponds to the bulk system free energy. For 
disordered systems, Eq. © is expected to hold, with the 
(configurationally-averaged) free energy given by Eq. 10} , 
and c taking the meaning of an effective conformal 
anomaly^. 

For our estimates of the effective central charge, we 
set T and p corresponding to the CNP, and took aver- 
ages of the free energy f(T c ,L) over three independent 
realizations. Figure shows the free energy at the CNP 
for 4 < L < 16, against 1/L 2 . A linear least-squares 
fit of the data gives c = 0.478(4). This is close to, but 
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FIG. 1: (Color online) Negative free energy at the conjectured 
Nishimori point, for strip widths 4 < L < 16 and N = 10 6 , 
against 1/L 2 . The line uses the central estimates from a 
parabolic least-squares fit to data, via Eq. 10, which gives 
f L = 1.7228(1) + 0.46(1)tt/6L 2 + 0.15(9)/L 4 . 

some 3% off, the result given in Ref. Q c = 0.464(4), 
which presumably was taken at those authors' own esti- 
mate of the location of the NP, p = 0.8906(2). On the 
other hand, the above result is compatible with the value 
corresponding to percolation in the Ising model, namely 
o p = 5V3 In 2/4tt w 0.4777i&. 

Incorporating curvature via the L~ 4 correction, as sug- 
gested by Eq. ©, results in the same f(T c , oo) (to within 
0.01%) as in the linear extrapolation. However, the con- 
formal anomaly estimate is changed to c = 0.46(1), which 
encompasses the result quoted in Ref. but appears in- 
compatible with Ising-model percolation. Earlier work in 
pureii and unfrustrated random-bondi£*A& Ising systems 
indicates that, for a square lattice, the L~ A term provides 
a important contribution towards stability and accuracy 
of free-energy extrapolations (note, however, that here 
the curvature effect is imperceptible to the naked eye, 
see Figure Should the same trend hold in the present 
case, c = 0.46(1) would appear to be the most reliable of 
the two estimates produced here. 



III. THE EXPONENT 77 

As a consequence of the preservation of conformal in- 
variance at a second-order phase transition, the correla- 
tion length £l on a strip geometry with periodic bound- 
ary conditions across (calculated at the critical point of 
the corresponding bulk two-dimensional system) is con- 



nected to the decay-of-correlations exponent 77, by the 
relationship^ 



For strips of homogeneous spin systems, the inverse of 
the dominant correlation length (related to the slowest- 
decaying critical correlations) is given by the first spectral 
gap of the transfer matria^S. A straightforward adap- 
tation for disordered cases can be devised through the 
replacement of pure-system eigenvalues by their counter- 
parts in a disordered environment, namely the Lyapunov 
characteristic exponentsiiiS^. One can then calculate 
the correlation length £z, (and thus the exponent 77 from 
Eq. (0), via 

CE 1 -aP-a« (8) 

where Af\ are, respectively, the largest and second- 
largest Lyapunov exponents. 

In the present case, the dominant correlations are 
ferromagnetic, and the Hamiltonian is invariant under 
global spin inversion. Therefore, in order to calculate 
A^ (A^ 1 '), it is sufficient to iterate \vq) (|i>i)) which is 
even (odd) under that same symmetrji2Si21, with no need 
for decontamination of the iterates of \v\). 

Before going further, one must recall that correlation 
functions at the NP are multifracta l 2 i 8 i 9 i 12 i 22 i 23 . In other 
words, the rate of decay (against distance R) of the 
moments of assorted orders, G n (R), of the correlation- 
function distribution is regulated by a set of exponents 
{i]n}, which are not connected by a single gap exponent, 
as is the case for pure systems where r) n = n r?i . The ex- 
ponent estimated via Eqs. Q and JHJ is in fact 770, which 
characterizes the zeroth-order moment of the correlation- 
function distribution, i.e. it gives the typical, or most 
probable, value of this quantity (see, e.g., Ref. and 
references therein). One has, in the bulk, 

G (R) = cxp [ln(a a R )] av ~ iT"° . (9) 

By evaluating estimates of 770 for the range of strip 
widths within reach of our computational facilities, we 
can, in principle, extrapolate the sequence L/tt^l to 
L — > 00, this way presumably accounting for higher-order 
finite-size corrections to Eq. 0. Earlier results for pure 
systems 1 ^ again indicate that L~ 2 is a convenient vari- 
able against which to set up an extrapolation scheme. 

We can also calculate corr elation functions directly on 
a strip, as done in Refs. I2l9lll2t and examine the behav- 
ior of their zeroth-order moment against distance, from 
which the appropriate correlation length can be extracted 
and plugged back into Eq. . Note that negative values 
of the correlation function will be present upon sampling; 
this is not an unsurmountable obstacle for the calculation 
of logarithmic averages here, as it is known that the dis- 
tribution at the NP is sharply peaked close to unitj£. 
Consequently, one can deal instead with absolute values, 
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FIG. 2: (Color online) Exponent T]l against 1/L 2 . Data are 
for L = 6, ■ • • , 16. The arrow pointing to the vertical axis 
indicates the percolation value, r/ p = 5/24. The straight line 
is a least-squares fit of L = 8 — 16 data. Shaded area close 
to the vertical axis indicates rough limits of confidence for 
extrapolation of finite-size data. Vertical bar to the left of 
vertical axis gives range of r/o from direct evaluation of zeroth 
moment of correlation-function distribution (see text). Verti- 
cal lines at left of graph show ranges of some recent estimates 
of rji. Full line, Refs. yj, \Q; short dashes, Ref. |l2j; long dashes, 
Ref.0 



as can seen by recalling that a logarithmic average is the 
same as the logarithm of a geometric mean: as long as the 
overall sign of the product of all terms is positive (which 
we have reason to believe here), it does not matter that 
some (few) are negative. 

In Figure El we present t)l calculated via Eqs. @ 
and (8J, against 1/L 2 . On increasing L from 4 (not 
shown) to 8, there is a decreasing trend in 77^ which ap- 
pears to halt, and turn to a roughly L— independent be- 
havior, when values corresponding to larger L < 16 are 
considered. Given the circumstances, the safest course of 
action is (i) to assume that the approximately constant 
behavior will not change significantly for larger L outside 
our computational capability range, and (ii) to extrapo- 
late the data at hand in as simple a manner as possible, 
treating the results with a large dose of skepticism. 

From a linear least-squares fit of L = 8—16 data 
against L~ 2 , shown in Figure we get a central ex- 
trapolated estimate 770 ss 0.192. By considering the error 
bars associated to finite-size estimates, it appears that 
any value in the range 0.187 < 770 < 0.196 (the extent of 
the shaded area in the Figure) would be plausible. 

We also evaluated correlation functions directly, and 
estimated the zeroth moment of their probability distri- 



bution. Following Ref. |9j, we used L — 10, and strip 
length TV = 10 7 columns; we found that the best set of 
results was for correlations calculated along the strip, for 
distances 1 < x < 18; when plotted on a semi-logarithmic 
scale, our data show slight curvature for 1 < x < 3, and 
set in to a very good straight line for 4 < x < 18, from 
which one gets t] = 0.194(1) (shown in the Figure, as 
a thick bar immediately to the left of the vertical axis) 
via Eq. . This is consistent with, and more accurate 
than, the extrapolation of the Lyapunov-exponent data 
given above. 

Comparison against data from previous work is as fol- 
lows. Numerical estimates from direct evaluation of the 
first moment of the probability distribution of spin-spin 
correlation functions give 771 = 0.182(5) (square lattice, 
approximate location of the NP at p = 0.8905(5)) 7 ; 
rji = 0.1854(19) (square lattice, approximate location of 
the NP at p = 0.8906(2))2,; ^ = 0.1854(17) (square lat- 
tice, CNP)&; 771 = 0.181(1) (triangular and honeycomb 
lattices)^. All are displayed in Figure for ease of vi- 
sualization. While overlap between these and the error 
bars of the present result is not better than marginal, 
it is clear that all estimates, for 770 and 771, exclude the 
percolation value, r\ v — 5/24 = 0.208333 • • • 25 (shown in 
the Figure by an arrow) by a safe gap. 

We sum up the situation as follows. Using Eq. (jHJ 
as a definition of the correlation length for random sys- 
tems is well justified in theoryi^G*^, and gives the in- 
verse decay rate of the zeroth moment of the correlation- 
function distribution. As seen above, in the present case 
the associated exponent 770 appears to differ slightly from 
771 which relates to the first moment. This is probably 
the rule rather than the exception; indeed, it has been 
shown that, for unfrustrated Ising systems, the finite- 
size scaling of numerical estimates of 770 derived in the 
context of Eqs. and © differs from that of results 
obtained directly from the spatial decay of correlation 
functionalist. Though in that case the origin of the dis- 
crepancy was traced to effects of the marginal disorder 
operator^ known to arise in the absence of frustration, 
the analogous operator structure at the NP is not known 
so far. However, it seems plausible to ascribe the small 
differences between the same two groups of results here, 
to similar causes. 



IV. SUSCEPTIBILITIES 

The uniform zero-field magnetic susceptibility \l on a 
strip is given by the second derivative of the free energy, 
relative to a uniform field h: 



XL 



d 2 f L 



dh 2 



(10) 



h=0 



As usual in the numerical calculation of derivatives, care 
must be taken to avoid introduction of spurious errors. 
We have considered an infinitesimal field 5h = 10~ 4 (in 
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units of Jo), for the finite differences used in the differ- 
entiation denoted in Eq. ijlOjl . We have also used the 
same configuration of bonds (that is, the same sequence 
of pseudorandom numbers) for the comparison of free en- 
ergies at different field values: free energies of the same 
bond geometry have to be subtracted. Thus, fluctua- 
tions in the finite differences used in the calculation of 
derivatives are much smaller than those for the free en- 
ergies themselves 7 *. For the calculations reported in this 
Section, we typically used strip lengths N — 10 7 . 

Finite-size scaling arguments suggest the following be- 
havior for xl at the critical point T c : 



XL 



(11) 



where 7 and v are, respectively, the exponents character- 
izing the singularities of bulk uniform susceptibility and 
correlation length. 

Another quantity of interest is the non-linear suscep- 
tibility Xz? 1 ^ ^ given in terms of the power-law expansion 
of the magnetization m: 



= xh-x (nl) h 3 



where 



X 



(nl) = 



dh 3 



d 4 P 



dh 4 



(12) 



(13) 



h=0 



The nonlinear susceptibility at criticality obeys a finite- 
size scaling relationship similar to Eq. itTTll . with the 
replacement 7 — > j nl . This quantity has been in- 
vestigated in the context of critical phenomena in both 
pur o 2 fi< 27 and (quantum) spin-glass magnets^. The nu- 
merical procedures described above, for the calculation 
of derivatives, are followed here as well. 

In Figure we show data for the linear susceptibility 
Xl, evaluated at the CNP, which have been fitted to the 
single power-law form, Eq. Ill 111 . We noticed that the 
X 2 per degree of freedom decreases from 2.3, for a fit 
including L = 4 data (from which the estimate j/v = 
1.82(1) is extracted), to 0.38 for a fit of L = 6 — 16 data 
only, thus pointing to a clear improvement in the quality 
of fit. The latter procedure gives j/u= 1.797(5), which, 
for the reasons just mentioned, we assume to be the best 
result to be extracted from the present data set. This is 
just consistent, at the margin, with the percolation value 
( 7 /i/) p = 43/24 « 1.7917^. 

In previous work, the following results have been 
found: j/v = 1.80(2)2 (square lattice, approximate lo- 
cation of the NP at p = 0.8905(5)), 7/1/ = 1.795(20)12 
(triangular lattice) . Both are consistent with the present 
estimate. 

Non- linear zero field susceptibility data (Xl), evalu- 
ated at the CNP, are exhibited in Figure A fitting 
procedure, similar to that used to extract the linear sus- 
ceptibility exponent, leads to 7"' jv = 5.59(2) when data 
for L = 4 — 16 are used, and to y nl jv = 5.55(2) when 
L = 4 data are discarded. However, the x 2 P er degree of 
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FIG. 3: (Color online) Double-logarithmic plot of zero-field 
linear susceptibility \l at the conjectured Nishimori point, 
for L = 4 — 16. The line is a single-power law least-squares fit 
to L = 6 — 16 data, enabling the estimation of y/v with use 
of Eq. (JllJ- 



freedom decreases only from 1.3 to 0.75 between the for- 
mer and latter fits. Bearing in mind that we are dealing 
with a small number of finite-size estimates, such a varia- 
tion does not warrant discarding L = 4 data on grounds 
of a significant improvement in the quality of fit. The 
value 7™/i/ = 5.59(2) is consistent with the scaling re- 
lation 7 ni jv = 2y/v + d = 5.59(1) (using yjv obtained 
above, and d — 2 for the space dimensionality). 



V. CONCLUSIONS 

We have calculated assorted critical quantities at the 
conjectured exact location of the Nishimori point (CNP) 
for square-lattice ±J Ising spin glasses, namely p — 
0.889972 • • • , T/Jq = 0.956729 • • • . By working at this 
fixed location, we attempt to eliminate one among many 
sources of uncertainty with which one has to deal in the 
study of disordered systems. Of course, whether or not 
such choice in fact introduces systematic distortions can 
only be found by comparison of a body of results per- 
taining to the problem under scrutiny. 

Our extrapolation of finite-size free-energy data, in 
order to produce an estimate of the effective central 
charge, has been careful by accounting for curvature ef- 
fects which are known to be relevant in such circum- 
stance a 1 ^ 17 ' 1 ^ . Our final estimate, c = 0.46(1), is con- 
sistent with an earlier result 2 , c = 0.464(4), calculated 
at p = 0.8906(2), and appears to exclude the percolation 
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FIG. 4: (Color online) Double-logarithmic plot of zero-field 
nonlinear susceptibility Xl 1 a t the conjectured Nishimori 
point. The curve is a single-power law least-squares fit to 
L = 4 — 16 data, enabling the estimation of 7 ni /^ with use of 
Eq. JIIJ- 

valued, c p « 0.4777. Had we not included curvature ef- 
fects, we would have reached c = 0.478(4) which would 
lead to the opposite conclusion. 

We conclude that whatever differences may exist be- 
tween free energies evaluated at the CNP, and those cal- 
culated at nearby locations such as that given in Ref. 0, 
their effect upon subsequent estimates of the central 
charge is not detectable amidst the noise associated to 
other sources of uncertainty. Prominent among these is 
the current upper limit on strip widths, L sa 20, imposed 
by practical considerations. 

We have evaluated finite-size correlation lengths via 
the difference between the two largest Lyapunov expo- 
nents. With the help of conformal-invariance concepts, 
these were used to produce a sequence of estimates of the 
decay-of-correlations exponent 770 , related to the decay of 
the zeroth moment of the correlation-function probabil- 
ity distribution. Though such a sequence does not behave 
as smoothly as its free-energy counterpart, it seems safe 
to state that it points to 0.187 < 77 < 0.196. We have 
also directly calculated correlation functions on a strip, 
thus assessing the statistics of the above-mentioned ze- 
roth moment. The corresponding result, 770 = 0.194(1) 
is consistent with, and more accurate than, that derived 
from the Lyapunov exponents. Both estimates slightly 
differ from 771, related to the first moment of the same 
distribution, for which available estimate o 2 i 7 i 9 i 12 fall in 
the range 0.180 < 77 < 0.187. In all cases, for 770 as well 
as 771, the random-percolation valued rj p = 0.208333 • • • , 



is definitely excluded. 

From zero-field susceptibility data we obtain 7/7/ = 
1.797(5), which falls within the range of previous re- 
sults 7 *^, and just about touches the percolation value 
(7/V) p = 43/24 w 1.7917-2&, at the lower end of the error 
bar. Similarly to the case discussed in Ref. fl^L it appears 
that from 7/1/ alone it is hard to get conclusive evidence, 
either for or against the behavior at the NP being in the 
percolation universality class. 

As regards non-linear susceptibilities, our study has 
been motivated by the well-known manifestations of 
multiscaling behavior of correlation functions at the 
]\jp2 i 8 i 9 i i2 i 22 i 23 rp^g connec tion between linear suscepti- 
bility x an d the first moment of the correlation-function 
distribution is given through the fluctuation-dissipation 
theorem, which (upon invoking standard scaling argu- 
ments 2 ^) implies the scaling relation 7/2/ = 2 — 771. The 
non-linear susceptibility X \ on the other hand, can be 
expressed in terms of four-point correlations and prod- 
ucts of two-point ones^. Thus it is not obvious a priori 
whether any of the multiscaling properties, observed for 
the assorted moments of the two-point function, will in- 
fluence x^ ■ Whatever guidance we have on the subject 
is given by the standard finite-size scaling relation be- 
tween the exponents associated to x an d X ■> namely 
-y"'/V = 27/V+d. This is established upon consideration 
of finite-size scaling properties of the free energy 2 ^, there- 
fore bypassing any explicit connection to correlation func- 
tions. Should multiscaling behavior of magnetization- like 
quantities occur (via their connections to aggregated cor- 
relation functions), one would expect to see something 
similar to the non-constant gap exponents observed for 
correlation function statistics (i.e., non-constant magne- 
tization gap exponents^) , which would imply breakdown 
of the relationship just mentioned. As seen above, we 
have found that the relationship is in fact obeyed, to 
very good numerical accuracy. Thus, no evidence has 
been detected for multiscaling behavior of magnetic sus- 
ceptibilities. 

NOTE ADDED. After initial submission of this paper, 
new work came up2& in which the results of Ref. Q are 
extended and reanalysed. The estimates of the location 
of the NP, and of the central charge, remain unchanged 
at p = 0.8906(2) and c = 0.464(4), respectively. 
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